Renormalized Jellium model for charge-stabilized colloidal suspensions 
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We introduce a renormalized Jellium model to calculate the equation of state for charged col- 
loidal suspensions. An almost perfect agreement with Monte Carlo simulations is found. Our 
self-consistent approach naturally allows to define the effective charge of particles at finite col- 
loidal density. Although this quantity may differ significantly from its counterpart obtained from 
the standard Poisson-Boltzmann cell approach, the osmotic pressures for both models are in good 
agreement. We argue that by construction, the effective charge obtained using the Jellium approx- 
imation is more appropriate to the study of colloidal interactions. We also discuss a possibility of a 
fluid-fluid critical point and show how the new equation of state can be used to shed light on the 
surprising results found in recent sedimentation experiments. 
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In spite of the great effort invested in trying to un- 
derstand the phase stability of colloidal suspensions, our 
knowledge of these complex systems is still quite rudi- 
mentary. It is curious to compare this situation with an 
earlier debate concerning the nature, or even the pos- 
sibility of the liquid-gas phase separation in symmetric 
electrolytes. Now this debate is almost over, and the 
phase structure of a symmetric 1:1 electrolyte is well elu- 
cidated, although the universality class of the critical 
point is still being discussed. The Coulombic critical- 
ity — to distinguish it from the solvophobic criticality — 
cannot be observed in water and organic solvents of low 
dielectric permittivity must be used 0. The phase sep- 
aration results from the strong electrostatic correlations 
between the cations and anions of electrolyte . 

The phase stability of charged colloidal suspensions, 
on the other hand, is far from being well understood 
Uti 1^ 1^ lZk >^ lEi JH] ■ The strong charge and size asym- 
metry between the macroions and the microions present 
in a suspension makes it very difficult to apply to these 
systems the traditional tools of the liquid state theory. 
Nevertheless, in view of results on symmetric 1:1 elec- 
trolytes, it would not be very surprising if colloidal sus- 
pensions inside a solvent of sufhcicntly low dielectric per- 
mittivity e also presented a gas-liquid phase transition. 
Indeed such an instability has been observed in recent 
Monte Carlo simulations, see for example and ref- 

erences therein. What is much more surprising is that 
there are some experimental indications of an instabil- 
ity even in aqueous suspensions containing only mono- 
valent counterions. Theoretical estimates of the strength 
of electrostatic correlations, for aqueous suspensions sug- 
gest that they should be too "hot" for an instability to 
set in. Nevertheless the experimental situation remains 
unclear [T^ . 

A number of theories have been proposed to address 
this unsettling experimental situation. A major draw- 
back of most of these approaches is that they rely on 
uncontrolled approximations which have not been fully 
tested. However, numerical solution of the full non- 



linear Poisson-Boltzmann (PB) equation inside a spheri- 
cal Wigner-Seitz (WS) cell finds no indication of any 
thermodynamic instability '7*1 . Of course, one can rightly 
question the reliability of the Wigner-Seitz cell PB model 
for the study of a fluid phase of a highly disordered sus- 
pension. On the one hand, since the main contribution to 
the osmotic pressure inside an aqueous suspension comes 
from the polyion-microions interactions, the cell might 
not be such a bad approximation. On the other hand, the 
cell fails to properly account for the colloid-colloid cor- 
relation, but these might not be of much importance for 
aqueous suspensions with monovalent counterions. Nev- 
ertheless, while the cell model is a good approximation 
for dense colloids, one should be very careful in extrap- 
olating its findings to highly dilute suspensions. Clearly 
there is an urgent need for an accurate theory which will 
not rely on the cell approximation and which would be 
relevant for the study of colloidal phase stability. In this 
Letter we shall present such a theory. Our approach is 
similar to the Jellium approximation much used in the 
solid state physics. 

Consider an aqueous suspension of colloidal particles 
of charge —Zq and radius a in contact with a reservoir 
of monovalent salt at concentration Cs and electrostatic 
potential 4>r = Q {q is the elementary charge). The num- 
ber of counterions and colons inside the suspension is de- 
termined by the thermodynamic equilibrium. While the 
colloidal particles are more or less uniformly distributed 
throughout the solution — we are mainly interested in 
the small density regime — , the positions of counteri- 
ons and colons are strongly correlated with the positions 
of polyions. As a leading order approximation we can, 
therefore, take the polyion-polyion correlation function 
to be Qpp = 1 |14|| while the exact polyion-counterion and 
polyion-coion correlation functions are gp± = g-P'^'±('>') ^ 
where w± (r) are the polyion-counterion and the polyion- 
coion potentials of mean force and /3 — \/{kBT) is the 
inverse temperature. 

Choosing the coordinate system in such a way that it 
is centered on top of one of the colloidal particles, the 
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electrostatic potential satisfies the Poisson equation 
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The charge density is pq{r) = -pback + qp+ir) - qp-{r), 
where, pback = Z]^a.ckqPp, and pp is the mean density of 
colloids inside suspension. The background charge is ex- 
cluded from the colloidal interior. Naively one can sup- 
pose that Zback = Z. This, however, is not correct and 
the bare charge must be renormalizcd in such a way as to 
lead to a self-consistent solution of Eq. as discussed 
below. 

Approximating the potential of mean force by the elec- 
trostatic potential, the local concentration of counterions 
and colons inside a suspension is, 
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It is important to keep in mind that in order for sus- 
pension to be neutral, the electrostatic potential at in- 
finity (bulk) cannot vanish but must saturate to a value 
(/)(oo) = 4>D given by Cs sinh[/3g(/)£,] = Pp.^back- There 
exists, therefore, an electrostatic potential difference be- 
tween the suspension and the salt reservoir. In the bio- 
physics literature, this potential difference is known as 
the Donnan potential, and is partially responsible for the 
biological cell trans-membrane potential. 

Far away from the colloidal surface, the electrostatic 
potential reduces to the familiar Debye-Hiickel expression 
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where — AnXBip+ioo) + p^{oo)] and As = Pq^ /f- de- 
notes the Bjerrum length. The value of Z^ti is determined 
self-consistently from the numerical solution of Eq. 
so that Zes{Z, Zhack,Cs) — Zhack- This renormalization 
of background charge is a consequence of counterion con- 
densation. It is important to keep in mind that for highly 
charged colloidal particles, ZeS is not equal to the bare 
colloidal charge. Furthermore, using the contact theorem 
and the vanishing of the electric field as r — > oo, the os- 
motic pressure within the suspension takes a simple form 



pp = pp+p+{oc)+p^ioo)-2c, - Pp+y/^>p + 4c2-2c, 

(4) 

It is a nice feature of the Jellium approximation that once 
the effective charge is determined, the Debye length 1/k 
and the osmotic pressure both follow directly. In Fig. 
^we compare the osmotic pressures calculated using the 
renormalizcd Jellium approximation to the ones obtained 
within the WS cell (for both approaches, the effective 
charges and pressures become independent of Z when the 
latter is large enough, approximately Z > 20a/ Xb', Fig. 
n has been obtained under this condition of saturation, 
which is usually met in colloidal suspensions). A surpris- 
ingly good agreement is found between the two theories, 




FIG. 1: Pressure- volume fraction isotherms within the 
Poisson-Boltzmann cell model (PBC) and the renormalized 
Jellium approximation for Cs = (deionized situation). The 
packing fraction is defined as = Aivppa^ /3. The inset shows 
the same data on the linear scale. 



with a discrepancy only for volume fractions 77 > 0.15. In 
Fig 121 the osmotic pressure calculated in the renormal- 
ized Jellium model is compared to the results of recent 
Monte Carlo simulations :4j (where a model system of 
charged spherical macroions and point counterions inter- 
acting solely through hard sphere and Coulomb forces 
has been considered). The agreement is excellent, and 
justifies the neglect of colloid-colloid correlations in our 
treatment. Clearly, for high packing fractions (namely 
77 > 0.1) colloid-colloid correlations become important 
and should invalidate our approach. 
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FIG. 2: Comparison between the osmotic coefficient 
PP/iZpp) = Zeff/Z calculated using Monte Carlo (MC) sim- 
ulations 3 and the renormalized Jellium model for Cs — 0, 
Z — 40, 77 = 0.00125. At large couplings, we predict 
(3P/{Zpp) 0.1802 a/As, shown as the dashed curve. 
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In Fig. the effective charges calculated using the 
cell model and the renormalized JeUium theory are 
compared. At very low volume fractions (r/ < fO^^) the 
saturated effective charge of a salt-free suspension is in 
perfect agreement with the value of saturated effective 
charge found within the WS cell and can be approxi- 
mated by a simple equation, 

^sat«T^[<5-7M??)] , (5) 

As 

where 7 « 1 and (5 « 2. For volume fractions i] > 10^''', 
there is a fairly strong disagreement between the effective 
charges predicted by the two theories, even if both models 
are in very good agreement for the value of the osmotic 
pressure. 
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FIG. 3: Effective values of colloidal charges calculated with 
the Poisson-Boltzmann cell model and the renormalized 
JeUium approximation for Ca = 0. Inset shows the same data 
but on a linear scale. 

If the effective charge is to be used to study the struc- 
tural properties of colloidal suspensions (such as the 
structure factors or any other quantity requiring the 
knowledge of the effective interaction potential between 
two polyions) , we argue that the JeUium effective charge 
is the more relevant quantity. The reason for this is that 
within the JeUium approximation, two colloids at large 
separations R interact by the usual DLVO potential, as 
follows after some algebra from integrating the stress ten- 
sor over a colloid's surface (see also 0) 

V{R) = gffg' . (6) 

e(I + na)"^ R 

This is not the case for two colloids inside a WS cell 
for which the interaction potential is a very complicated 
function of separation and can only be calculated nu- 
merically. In fact, for separations larger than the inter- 
particle distance, it is even difficult to properly define 
what one means by an interaction potential between col- 
loids within the WS formalism, since the overall charge 



neutrality results in a vanishing electrostatic interaction 
between the two cells. 

To assess the quantitative validity of our approach, we 
compare our (saturated) effective charges to those de- 
duced in recent experiments using a torsional resonance 
spectroscopy 0. For colloidal volume fractions in the 
interval 10~^ < r] < 3.10~^, ZgatAs/a was found to be 
close to 6 This is in much better agreement with 

the renormalized JeUium model, which finds for these 
volume fractions ZgatAs/a varying between 6.7 — 7.0, as 
compared to the WS cell prediction of 8.7 — 9. We also 
note that a value close to 6 has been reported for this 
density regime in the Monte Carlo study presented in 
Ref [13. 

We also emphasize that within the JeUium model, the 
inverse screening length k is naturally related to the effec- 
tive salt density [see the definition below Eq. ©]. This 
should be contrasted with the cell model, for which there 
is no simple connection between the two quantities. 

The JeUium model predicts stability of a charged col- 
loidal suspension against a fluid-fluid phase separation. 
This reinforces the cell picture, where no instability is 
found. Our approach, however, neglects microionic cor- 
relations, which become important at high electrostatic 
couplings 2, 18., .19] ■ more precisely when F = v^Xs/d 
exceeds a threshold close to 2, d being a characteristic 
distance between microions in the electric double-layer 
and V their valency. In practice, however, ionic hydra- 
tion puts a lower bound to d, which prevents the high 
coupling regime from ever being reached in water with 
monovalent counterions. Under these conditions our ap- 
proach should, therefore, be quite reliable. Alternatively, 
when dealing with point particles, d may be estimated as 
d ~ {Ana'^v/Zy/'^ [lliSj- For the Monte Carlo simula- 
tions reported in Fig. [3 F = 2 therefore corresponds 
to an instability threshold As /a ~ 1.1. Beyond this 
point, the microionic correlations destabilize the system 
and lead to a fluid-fluid separation 0, 0|. Neverthe- 
less, the pressures obtained within the renormalized Jel- 
lium approximation are in very good agreement with the 
Monte Carlo simulations, even at fairly large electrostatic 
couplings (in Fig. [3 we have included the highest cou- 
pling for which the pressure was computed in Q). We 
therefore expect that the JeUium approximation, suit- 
ably corrected to include the counterion induced attrac- 
tion between the colloidal particles present at strong cou- 
plings, might be sufficient to account for the thermody- 
namic instability in suspensions containing multivalent 
ions. Work along these lines is in progress. 

Finally, it is interesting to speculate how the renor- 
malized JeUium model can help to shed new l igh t on the 
problem of sedimenting colloidal dispersions "2^, "211,121 
lil J4]. Recently, Philipse and Kocnderink (PK) 25j ob- 
served strongly inflated concentration profiles for charged 
mono-dispersed colloidal particles in absolute ethanol. 
Use of ethanol, instead of water, allows to produce highly 
deionized suspensions with "salt" concentrations on the 
order of 10"^ M. The renormalized JeUium predicts that 



4 



for this salt concentration, an infinitely dilute colloid 
would have an effective charge of Z^g « 20 a/ As . It 
is difficult to know what is precisely the bare charge of 
colloidal particles inside ethanol, however since the Bjer- 
rum length in ethanol is Xb ~ 2.3 nm, a fairly small bare 
charge on the order of a few hundred electrons is enough 
to place the particles (with radius a ~ 90 nm 25]) in the 
saturation regime. Using this observation we can par- 
tially account for the observations of PK. 

Static equilibrium of a suspension in a gravitational 
field requires that 

^ = -^9Pr> , (7) 

where m is the colloidal mass (corrected for buoyancy), 
g the gravitational acceleration, and x the vertical dis- 
placement. It is convenient to define the gravitational 
length as Ig = l/{mqj3), which for the experiments of 
PK is Ig « 0.2 mm j25i]. Substituting the expression for 
the osmotic pressure Eq.Q into Eq.Q and using the 
dependence of the saturated effective charge on the col- 
loidal volume fraction, Eq.© can be integrated. We find 
that for colloidal volume fractions rj > rji ~ 10^^ (i.e. in 
the counterion dominated regime): 

Mv/Vo) - ^ [ln^(»/) - ln^(77o)] = , (8) 

where r]o is the reference volume fraction at the x ori- 
gin. This can be taken as rjo « 0.01, the point at 
which the hard-core effects are completely negligible, and 
one can be sure to be looking at the dilute tail of the 
sedimenting profile. Note that for small volume frac- 
tions, rj ^ cxp{~ y/2XBx/alg) unlike the simple expo- 
nential barometric law. For rj < -qx « 5 x 10~®, the 
salt resulting from solvent dissociation dominates over 



counterions, and we recover the usual barometric law 
ri(x) ~ exp(— x/^g). Results concerning the crossover 
regime rjx < t] < rji will be published elsewhere. It is in- 
teresting to estimate the extent of the sedimentation pro- 
file predicted by the renormalized Jellium model. Since 
the concentrations of colloids on the order of rjm « 10~^ 
are still detectable, we find that the profile extends dis- 
tance X '-^ alg\n^ {T]m) / {2\b) « 2600 /g. For the colloids 
used by PK, this is almost 50 cm! Hence, the observed 
inflation of the sedimentation profile. 

To conclude, we have proposed a self-consistent renor- 
malized non-linear Jellium model that constitutes an al- 
ternative to the widely used PB cell approach for charged 
colloidal suspensions. Surprisingly, the equations of 
state within the two theories turn out to be very close, 
even though that the approximations involved to account 
for finite colloidal density are very different (finite cell 
against a renormalized background). Our results point 
to the relevance of the cell picture even at extremely low 
densities. This reinforces the argument that a fiuid-fluid 
phase instability is impossible, for highly charged colloids 
in water at room temperature, as long as the microions 
are monovalent [which allows to neglect microionic cor- 
relations, thereby identifying the potential of mean force 
with the electrostatic potential, see Eq. 0]. Our ap- 
proach also allows to define in a natural way, not only 
the effective charge of the macroions, but also the ef- 
fective screening length; these quantities are experimen- 
tally and conceptually more relevant than those obtained 
within the cell approach. Finally, our equation of state 
for dcionized systems helps in understanding recent sed- 
imentation experiments where "anomalous" density pro- 
files have been reported. 
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